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ABSTRACT 

Aims. We analyze the robustness of H-deficient post-AGB tracks regarding previous evolution of their progenitor stars 
and the constitutive physics of the remnants. Our motivation is a recent suggestion of Werner & Herwig (2006) that 
previous evolution should be important in shaping the final post-AGB track and the persisting discrepancy between 
asteroseismological and spectroscopical mass determinations. This work is thus complementary to our previous work 
(Miller Bertolami & Althaus 2006) and intends to shed some light on the uncertainty behind the evolutionary tracks 
presented there. 

Methods. We compute full evolutionary models for PG1159 stars taking into account different extramixing (overshooting) 
efficiencies and lifetimes on the TP-AGB during the progenitor evolution. We also assess the effect of possible differences 
in the opacities and equation of state by artificially changing them before the PG1159 stage. Also comparisons are made 
with the few H-deficient post-AGB tracks available in the literature. 

Results. Contrary to our expectations, we found that previous evolution is not a main factor in shaping H-deficient 
post-AGB tracks. Interestingly enough, we find that only an increase of ~ 50% in the intershell opacities at high 
effective temperatures may affect the tracks as to reconcile spectroscopic and asteroseismologic mass determinations. 
This forces us to conclude that our previous tracks (Miller Bertolami & Althaus 2006) are robust enough as to be 
used for spectroscopic mass determinations, unless opacities in the intershell region are substantially different. Our 
results, then, call for an analysis of possible systematics in the usually adopted asteroseismological mass determination 
methods. 
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, 1. Introduction 

. Post Asymptotic Giant Branch (AGB) stars constitute a 
' short-lived transition stage between AGB stars and white 
, dwarf stars. Among them a minority show H-deficient com- 
' positions and are suppossed to be the main progenitors of 
] H-deficient white dwarfs, which account for about 15% of 
the white dwarf population (Eisenstein et al. 2006). The 
group of H-deficient post-AGB stars displays a wide vari- 
ety of surface chemical compositions ranging from almost 
pure helium envelopes to the helium- (He), carbon- (C) and 
oxygen- (O) rich surface composition of the Wolf Rayet cen- 
tral stars of planetary nebulae ([WC]) and the PG1159 type 
stars; see Werner & Herwig (2006), from now on WH06. 
The surface composition of the last group resembles the in- 
tershell region chemistry of AGB star models when some 
overshooting at the base of the pulse driven convective zone 
(PDCZ) is allowed during the thermal pulses (Herwig et al. 
1997). For this reason, and also due to the fact that the oc- 
currence of late (i.e. post-AGB) thermal pulses is statisti- 
cally unavoidable in single stellar evolution modeling (Iben 
et al. 1983), a late helium shell flash is the most accepted 
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mechanism for the formation of these stars (see, however, 
De Marco 2002). In this scenario, the remaining thin H-rich 
envelope is either burnt in a very late helium flash (VLTP) 
that occurs on the hot white dwarf cooling branch after H 
burning has almost ceased, or diluted in a late helium flash 
(LTP) that develops when the H burning shell is still ac- 
tive during the horizontal evolution of the stars in the HR 
diagram (Herwig 2001). 

Roughly a third of PG1159-type stars exhibit multiperi- 
odic luminosity variations caused by non-radial g-mode 
pulsations. This has allowed researchers to derive struc- 
tural parameters — in particular the mass of these stars 
— of individual pulsators by means of asteroseismological 
studies i.e. by comparing adiabatic pulsation periods with 
the observed ones — e.g. Kawaler & Bradley (1994) and 
more recently Corsico & Althaus (2006). It is important to 
mention that for applications requiring accurate values of 
adiabatic pulsation periods full evolutionary models with a 
realistic thermal structure should be used. Stellar masses 
of PG1159 stars can also be derived by comparing the val- 
ues of log g and log Toff coming from the fitting of line- 
blanketed non-LTE model atmospheres to the measured 
spectra (Werner et al. 1991) with tracks coming from stellar 
evolution modeling. These two different approaches enable 
us to compare the derived stellar masses. Although previous 
spectroscopical mass determinations, based on old H~rich 
post-AGB models, show relatively good agreement with 
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astcroscismological masses (to about 5%, WH06, roughly 
0.03 Mq), the development of a new generation of stel- 
lar evolution sequences that account for the C- and O- 
rich surface abimdanccs expected in PG1159 stars (Hcrwig 
et al. 1999) has changed the situation. As mentioned by 
WH06 the new post-AGB tracks are systematically hot- 
ter than the old ones, which leads to lower spectroscopi- 
cal masses. The new mean spectroscopical mass becomes 
0.573 Mq, this is 0.044 Mq lower than previous values; see 
Miller Bertolami & Althaus (2006), from now on MA06. 
This is at variance with astcroscismological predictions. In 
fact from Table 3 of WH06 and Table 2 of MA06 the as- 
teroseismological masses are usually 10% higher than their 
spectroscopical counterparts, except for the hottest known 
pulsating PG1159 star RX J2117.1-I-3412, the spectroscop- 
ical mass of which is more than 20% higher than the as- 
tcroscismological one: sec Corsico ct al. (2007) for a recent 
and detailed study of this object. The difference in derived 
masses is a clear indication of the uncertainties weighting 
upon the mass determination methods. 

In this context, WII06 have recently compared new and 
old tracks and claimed that the previous evolution on the 
thermally pulsing AGB (TP-AGB) — particularly the third 
dredge-up (3DUP) efficiency and mass-loss rates — plays 
a decisive role in the location of the tracks in the HR and 
log Tcff— logg diagrams during the post AGB evolution. 
Specifically, as shown by Herwig et al. (1998), a strong 
3DUP changes the evolution of the core mass without al- 
tering the evolution of its radius. Consequently the mass- 
radius relation of the remnants will depend on the previ- 
ous TP-AGB evolution and, if we accept in the "predic- 
tion" of shell homology relations (igheii ~ -^core-^core""^, 
for Mcore ^ O.8M0, Herwig et al. 1998), the post-AGB 
tracks would be accordingly altered. WH06 also point out 
that mass loss can produce a similar effect as remnants of 
similar mass may come out with very different degrees of 
degeneracy depending on the previous evolution. This being 
the both mass-loss rates and 3DUP efficiency are 

poorly known, the location of theoretical post AGB tracks, 
and thus mass determination, would be highly model de- 
pendent and uncertain. These issues call for the need of 
an analysis of the robustness of existing H-deficient post- 
AGB tracks and for a way of solving the mentioned mass 
discrepancies 

However, no calculation of the importance of these ef- 
fects was actually presented neither in WH06 nor in Herwig 
et al. (2006). The lack of consistent calculations to assess to 
what a degree the location of the post AGB tracks depends 
on the prior AGB evolution has motivated us to undertake 
the present investigation. In the following sections we will 
elaborate on these issues. In this sense the present work 
is complementary to that of MA06 where H-deficient post- 
AGB tracks were presented but no analysis of its robustness 
was performed. In Sect. 2 we analyze how evolution pre- 
vious to the PG1159 stage affects PG1159 tracks in light 
of the suggestion presented by WH06. Then in Sect. 3 we 
explore to what an extent the constitutive physics of the 
models at the PG1159 stage may affect the tracks. In Sect. 
4 we compare with other H-deficient tracks available in the 
literature and also compare the location of tracks coming 
from LTP and VLTP events. Finally Sect. 5 is devoted to 
the discussion of the results and making some final remarks. 



2. Influence of previous evolution 

As was mentioned, uncertainties in mass-loss rates are ex- 
pected to affect the duration of the TP-AGB phase and to 
lead to remnants with different degrees of degeneracy and 
mass luminosity relation. Also the initial-final mass rela- 
tion is expected to be altered by different mass-loss rates. 
By altering the intensity of mass loss we can get the same 
final remnant mass from progenitors of initially very differ- 
ent mass, and consequently very different previous evolu- 
tions (e.g., that have or have not undergone a helium core 
flash at the tip of the RGB). We will address these issues in 
Sect. 2.1. In Sect. 2.2 we elaborate on the effects of differ- 
ent 3DUP efficiencies on the TP-AGB, which is the other 
point mentioned in WH06 as a possible cause for shifts 
in post-AGB tracks. The main effect of 3DUP efficiency 
is to change the initial-final mass relation. Indeed, strong 
dredge up events on the TP-AGB lead to lower final rem- 
nant masses for the same initial mass. In this context we 
analyze sequences with different 3DUP efficiencies and, to 
disentangle this effect from the one studied in Sect 2.1, with 
the same TP-AGB lifetime. Finally, Sect 2.3 is intended 
to clarify the reason of the difference between MA06 and 
Blocker (1995a, b) tracks and to study the extreme limiting 
case for which no overshooting (OV) is allowed to operate 
at any convective boundary during the whole evolution. In 
all the sequences presented in this section, mass loss has 
been arbitrarily set during the departure from the AGB 
as to get a VLTP and the subsequent PG1159-like surface 
composition. 

To visualize and quantify the change introduced by the 
variations in the parameters of each sequence we will refer 
and compare our sequences with those of MA06 and Corsico 
et al. (2006) which are assumed as standard in the present 
work. These sequences were calculated with an overshoot- 
ing efficiency of / = 0.016 at all convective boundaries; see 
Herwig et al. (1997) for a definition of /. To quantify the 
change that a variation in T^g and g for a sequence of a 
given mass — caused by different physical assumptions in 
the calculations would produce in spectroscopical mass 
determinations, we estimate a mass value for the sequence 
from its location relative to MA06 sequences this is what 
is called "mass derived from comparison" in Table 1 — and 
compare that mass with the actual value of the mass. The 
difference between both values gives the shift expected in 
spectroscopical masses if tracks with a different physical 
assumption are used in their derivations. 

It is worth noting that most of our article deals with 
post- VLTP sequences. However, within the late helium 
flash scenario for explaining the origin of PG1159 stars, 
these objects could also be the offspring of LTP events. 
In fact some PG1159 stars are known to be ^^N-deficient, 
a fact usually asociated with post-LTP objects. In these 
cases some H will be present but hidden below the detec- 
tion limit. If systematic differences exist between post-LTP 
and post- VLTP tracks, then this will introduce a system- 
atic effect in spectroscopic mass determiations. Although 
from figure 1 of Herwig one is tempted to conclude that 
this is not the case, it is worth noting that the presence of 
H should be more important in the low mass region as these 
stars display thicker H-envelopes. We will discuss this issue 
in section 4.3, where detailed comparisons between post- 
LTP and post- VLTP tracks will be made for a wide range 
of masses and various surface H-abundances. . 
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Sequence 


Final Mass 


Mass derived 
from comparison 


9 3 


- 
- 


NALT 


0.607 


0.612 


9.25 




T A T T 


U.UUOu 




- 


SALT 


0.6033 


0.598 






2.2MSALT 


0.5157 


0.524 


9.2 




TPA008 


0.617 


0.621 






TPA004 


0.633 


0.637 






3M0 w/NOV 


0.626 


0.623 


g 9.15 





Table 1. Values of the final masses of the sequences of this 
work and the masses derived from the comparison with the 
"standard" ones (MA06). Stellar masses are in M©. See 
Sect. 2.1, 2.2, and 2.3 for definition of the sequences. 



2.1. Effect of different TP-AGB lifetime or mass-loss rate 

As stated in WH06, for a similar core mass, a remnant 
that spend more time on the TP-AGB will finish with a 
more compact and degenerate core. Then, different mass- 
loss prescriptions can lead to remnants with the same core 
mass but different radius and, by virtue of shell homology 
relations — that "predict" ighoU ~ -^corc-^coro: Herwig et 
al. (1998) — different luminosities. This is supported by 
the work of Herwig et al. (1998) that shows that, because 
of the continuous shrinking of the H-free core (HFC), the 
luminosity at the TP-AGB keeps increasing, despite the 
end of the effective core mass growth as consequence of 
strong dredge up events. In addition. Blocker (1995b) has 
already shown that a more compact remnant is more lumi- 
nous than a less compact one of similar mass. To analyze 
the possible shift in the H-dcficient post-AGB tracks result- 
ing from uncertainties in TP-AGB mass loss — and hence 
in different TP-AGB lifetimes — we have calculated the full 
evolution of three sequences with the same prescriptions as 
in MA06 but changing mass loss at the TP-AGB to get dif- 
ferent TP-AGB lifetimes. These sequences are: NALT with 
a normal mass-loss prescription, SALT with a short TP- 
AGB lifetime and LALT with a longer TP-AGB. AU these 
sequences come from the same pre TP-AGB evolution of an 
initially 3-Mq ZAMS star. While NALT underwent 12 ther- 
mal pulses, SALT and LALT sequences experienced 6 and 
18 pulses, respectively. SALT (LALT) sequence has a TP- 
AGB hfetime a factor 2 shorter (1.5 longer) than NALT. 
Thus the sequences considered here take into account possi- 
ble uncertainties in TP-AGB lifetimes up to a factor three. 
This is more than what is expected from different mass- 
loss prescriptions (Kitsikis & Weiss 2007). Due to the high 
dredge up efficiency during the last thermal pulses the fi- 
nal remnant mass of all these sequences is very similar (see 
Table 1), thus allowing a direct comparison of the effect of 
different TP-AGB lifetimes on the location of H-dcficient 
post-AGB tracks of similar mass. We mention that all of 
these sequences have been followed through an additional 
post-AGB thermal pulse (the VLTP) where the H-rich en- 
velope is violently burned. 

In some agreement with Blocker (1995b) we find a shift 
in post-AGB tracks as a consequence of different TP-AGB 
lifetimes. However the effect is not very important. In fact, 
comparing SALT and LALT sequences (both with the same 
final mass) we see that a factor 3 in TP-AGB lifetimes leads 
to a maximum shift of 0.03 dex in log T^s- A shorter TP- 
AGB leads to cooler tracks that would imply ~ 0.015 Mq 
larger spectroscopical masses. It is also worth noting that 
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Fig. 1. Evolution of the HFC (mass and radius) during 
the TP-AGB (solid lines, filled circles) and at the VLTP 
(dashed lines, empty circles) for selected sequences (values 
are taken just before each thermal pulse). Also the locus of 
the standard models at the moment just before the VLTP 
is shown for comparison. Note that, due to the turn off of 
the H-burning shell, compression before the VLTP does not 
follow the trend in the AGB. [color figure only available in 
the electronic version! 



tracks for LALT and NALT sequences are almost identi- 
cal regardless the difference in TP-AGB lifetime of 50%. It 
seems that, while shortening the TP-AGB does change the 
post-VLTP tracks, prolonging it does not produce a size- 
able effect. To understand this, we show in Figs. [1] and [2] the 
mass-radius relations of our sequence for both the H- and 
He- free cores — HFC and HeFC, respectively. The evo- 
lution of the HFC is in agreement with that presented by 
Herwig et al. (1998) which shows that the HFC continues to 
contract even when the core mass growth is stopped by ef- 
ficient 3DUP events. Because this behaviour is the basis of 
the argument of WH06 the following should be noted. First 
the radius of the HFC at the moment of the VLTP does not 
follow the trend during the TP-AGB. This is a result of the 
accelerated compression of the intershell caused by the de- 
cline of the H-burning shell when the star approaches the 
white dwarf cooling track. Second, and more importantly, 
the post-VLTP sequences are powered by the He-burning 
shell and consequently, if shell homology are to be used 
in the analyzes, the relevant values should be the HeFC 
mass and radius. Note that the HeFC (Fig. [5]) seems to 
converge to a certain locus in the core mass-radius diagram 
faster than the HFC. In fact while in all the sequences the 
HFC radius gets smaller with each thermal pulse, the HeFC 
ends its compression after about ^ 10 thermal pulses. This 
helps to understand why there is almost no difference be- 
tween NALT and LALT sequences. The 6 "extra" thermal 
pulses of LALT sequence do not introduce any significant 
change in the mass-radius relation and thus, according to 
shell homology relations, their post-AGB luminosity should 
be similar. 

As mentioned early, different mass-loss rates can also 
change the initial-final mass relation of the sequences, lead- 
ing to final remnants with very different previous evolution 
but similar final mass. In this connection, we have com- 
puted the evolution of an initially 2.2-Mq sequence by as- 
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Fig. 2. Same as Fig. [l] but for the HeFC. Note that the 
location of the HeFC on this diagram seems to converge, 
after not many thermal pulses, to a certain locus, [color 
figure only available in the electronic version] 



suming an extreme mass-loss rate during the whole AGB 
(sequence 2.2MSALT). As a result, this sequence underwent 
only 5 thermal pulses on the AGB — as compared with the 
15 AGB pulses of the 2.2-M0 sequence in MA06. The final 
mass of the remnant is of 0.516 Mq, much lower than the 
0.565 Mq quoted in MA06. The track for this sequence in 
the log Teff — log g plane is shown in Fig. [Altogether with the 
other sequences of this work and those of MA06. Note that 
the 2.2MSALT track is more luminous and hotter than that 
of the standard sequence of similar mass (the 0.512 Mq se- 
quence in MA06). Note that in this case, the shift in the 
M-g-Tcs relation of the remnants would imply a decrease 
of ~ 0.01 Mq in spectroscopical masses. This value is un- 
expectedly low in view of the fact that the two standard 
sequences in the same region of the log Teg— log g diagram 
have a very different previous evolution. Indeed, the 0.512 
and 0.53 Mq sequences in MA06 have been calculated from 
an initially 1-Mq progenitor that went through the helium 
core flash. Also, the 0.512 Mq has a very different intershell 
and surface composition with only 2%, by mass, of oxygen; 
see MA06 for a description of this sequence. Again, it is in- 
teresting to look at the structure of the HeFC to understand 
this change. As can be seen in Fig. [2] (black star symbol), 
although the mass and radius of this model fall almost in 
the standard locus (the thick grey line in figure 2), its HeFC 
mass (^ 0.45 Mq) is significantly higher than that of the 
standard 0.512 Mq sequence 0.43 Mq) and thus should 
be, again from shell homology arguments, more luminous 
than the standard sequence. Indeed that is what actually 
happens. Even more, the 2.2MSALT sequence has a HeFC 
mass that falls almost in the middle of that of the 0.512 and 
0.53 Mq MA06 sequences and its track in the log Teg— log g 
diagram does exactly the same. These considerations seem 
to support the idea that is the HeFC structure — and not 
the HFC — which is important to understand H-deficient 
post" AGB tracks. 

So, although the findings of this section confirm that 
different TP-AGB lifetimes may result in changes in the 
post-AGB tracks, we find that this effect is not enough to 
account for the mass discrepancy mentioned in the intro- 
duction. Indeed, we find that the PG1159 spectroscopical 




Log(T„) 



Fig. 3. PG1159 tracks of this work as compared with 
those of MA06. Thin solid lines correspond to the stan- 
dard (/=0.016 at all convective borders) tracks of MA06 
with stellar masses of (from right to left) 0.512, 0.53, 0.542, 
0.565, 0.585, 0.609, 0.664 Mq. 



masses inferred from the MA06 post-AGB tracks would be 
higher by at most ~ 0.015 Mq (for stars close to the 0.6 
Mq tracks) if in their calculations MA06 had considered 
much shorter TP-AGB lifetimes during the progenitor evo- 
lution of their PG1159 sequences. On the other hand we 
find impossible to get a similar shift for stars close to the 
0.512 Mq track. This is so because the lack of 3DUP in low 
mass stars. 

2.2. Effect of the third dredge-up efficiency 

To explore the role of the 3DUP efficiency during the TP- 
AGB in the location of post-AGB tracks, we have followed 
the TP-AGB evolution for three different values of the over- 
shooting efficiency (/) at the pulse driven convection zone 
(PDCZ) that develops during each He-shell flash. As shown 
in Herwig (2000), higher / values at the bottom of the 
PDCZ lead to more intense helium shell flashes and more in- 
tense third dredge up events, while the value of / at the base 
of the convective envelope only plays a secondary role in de- 
termining the 3DUP efficiency (the reasons for this are ex- 
tensively discussed and shown in sections 4 and 5 of Herwig 
2000). We have, thus, calculated three different sequences 
for a Z-Mq progenitor by adopting values of /=0.016, 0.008, 
0.004 at both convective borders of the PDCZ, from now 
on sequences NALT, TPA008 and TPA004; sequence NALT 
corresponds to that previously described. At any other con- 
vective zone — for example the AGB convective envelope 
and the core burning regions in the previous evolution — 
the "standard" value of /=0.016 has been used. We stress 
that the "standard" value / = 0.016 comes from the fitting 
of the width of the main sequence (Herwig et al. 1997), 
and thus is appropriate for the core H-burning zone. But it 
may be unrealistic for the conditions at the PDCZ (Herwig 
2004). All of these sequences have similar TP-AGB life- 
times. This enables us to disentagle the 3DUP effect from 
the one coming from different TP-AGB lifetimes, which 
was studied in the previous section. Also, for these three 
sequences, mass loss during the last interpulse phase has 
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been artificially set in order to obtain a VLTP and conse- 
quently a H-deficient post-AGB remnant. 
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Fig. 4. HFC evolution during the TP-AGB for three se- 
quences with different / values at the PDCZ (masses in 
Mo). 

Note from Fig. 2] that different values of / yields differ- 
ent evolution of the HFC. For models with higher 3DUP 
efficiencies the "effective" growth of the HFC is stopped. 
This is because the increase in the HFC induced by the H- 
burning shell is compensated for by a decrease during the 
3DUP events. Not only the HFC mass is altered but also, 
as expected, the surface and intershell abundances — in 
particular the O intershell abundance; see Herwig (2000). 
As a result of the different adopted 3DUP efficiencies, the 
final remnant masses are different, being 0.607, 0.617 and 
0.633 Mq for NALT, TPA008 and TPA004 respectively. It 
is worth noting that TPA004 hardly undergoes any 3DUP 
events. So this sequence should be representative of the case 
in which no overshooting is considered at the PDCZ. 

Our results suggest that different 3DUP efficiencies do 
not seem to lead to an important shift in the location of 
the theoretical post-AGB models in the M-g-T^s space. 
Indeed, sequences TPA008 and TPA004 are located in the 
zone of the log Teg— log g diagram corresponding to rem- 
nants of similar mass of the standard sequences; see Fig. [31 
A quantitative measure of the possible shift in the tracks 
relative to the standard MA06 ones is given Table 1. Note 
that there is a small shift of 0.005 Af© in the derived mass 
for the NALT sequence as compared to the actual one — we 
remind that NALT sequence has the same overshooting pre- 
scription than that assumed in MA06. This is probably due 
to a combined effect of a different number of thermal pulses 
and slightly different envelope composition — which leads 
to different intershell opacities, see Sect. 3. Because the 
three sequences have similar TP-AGB lifetimes, this small 
shift should be taken as the level of uncertainties in these 
comparisons. Keeping this in mind, the masses derived for 
TPA008 and TPA004 are practically similar to the actual 
masses of these sequences. This shows that, at least around 
0.6 Mq, the theoretical M-g-T^s relation of the MA06 H- 
deficient post- VLTP sequences does not seem to depend on 
the intensity of 3DUP events. This can be understood from 
Fig. [2l Note that HeFC mass-radius values of sequences 
TPA008 and TPA004 at the moment of the VLTP lie on 



the same locus than the standard MA06 sequences of sim- 
ilar masses. Thus, according to shell homology relations, 
the He-shell luminosity-mass relation for these sequences 
should be similar to the MA06 ones — which do experience 
efficient 3DUP events. Finally, we mention that the central 
values of density and temperature {Tc, Pc) show that the 
HeFC readjusts its structure to the new mass after each 
thermal pulse. At the end of the TP-AGB the Tc, pc values 
of TPA004 are within those of NALT — of final HeFC mass 
0.572 Mq— and those of the 3.5 Mq sequence of MA06 — 
of final HeFC mass 0.638 Mq — , a fact which is consistent 
with the final HeFC mass of 0.601 Mq that characterizes 
sequence TPA004. 



2.3. Evolution previous to tlie AGB 



We explore now the effect of overshooting efficiency dur- 
ing both the early AGB and the core He-burning phase 
on the location of the post-AGB tracks. This bears also 
some relevance on the fact that, as inferred from the two 
previous sections, neither the TP-AGB lifetime nor the 
3DUP efficiency are the reasons for the fact that the MA06 
tracks are markedly hotter than the older H-rich tracks 
(Blocker 1995b). To assess these issues, we have calculated 
the evolution of an initially 3-Mq progenitor but without 
overshooting mixing at any convective border of the star 
during its whole evolution. After 19 thermal pulses, a H- 
deficient post-VLTP sequence of 0.626 Mq is obtained — 
early AGB and TP-AGB Hfetimes are ~ 4.8 x 10''yr and 
~ 2.1 X lO^yr, respectively. This is similar to one of the 
sequences of Blocker (1995b) that consisted of an initially 
3-A/o model that after 20 thermal pulses ends its post- 
AGB evolution as a 0.625 Mq remnant — with early AGB 
and TP-AGB lifetimes of - 7. x lO^yr and - 1.9 x lO^yr, 
respectively — and will allow us for comparison. The main 
evolutionary difference between both sequences is the oc- 
currence of a VLTP in the post-AGB evolution of our se- 
quence. 

The resulting H-deficient post-VLTP track is very sim- 
ilar to the MA06 one and thus much hotter than the old, 
H-rich, Blocker's O.626-M0 track. In fact if we estimate its 
mean mass from comparison with the standard MA06 se- 
quences we get almost the actual mass (see Table 1). The 
mayor difference is that the model is slightly cooler at the 
knee — a shift that would affect spectroscopical masses less 
than ~ 0.01 Mq. From Fig.[2]we can see that the evolution 
of the HeFC mass-radius relation is different from that of 
the standard sequences. But even in this case the differ- 
ence in the radius of the He- free core at the moment of the 
VLTP amounts to only a 4% as compared with the stan- 
dard sequences of similar mass. Consequently, it should not 
be surprising that the tracks are similar. 

This shows that the M-g-T^s relation for the post-AGB 
tracks is not significantly affected by the previous evolution. 
Thus, differences in the previous evolution do not seem to 
provide a possible solution to discrepancy between aster- 
oseismological and spectroscopical masses nor an explana- 
tion to the difference with older tracks. 
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3. The role of microphysics and composition in the 
location of post-VLTP tracks 

We explore now the importance of microphysics and chemi- 
cal compositions. Specifically we assess the effects of chang- 
ing the equation of state (EoS), chemical composition of 
the C-0 core and opacities — both radiative and conduc- 
tive. Here, we do not calculate new evolutionary sequences 
from the ZAMS to the PG1159 stage; instead we consider 
some post-VLTP sequences of MA06 and alter their micro- 
physics before entering the PG1159 stage. We have checked 
that the models are already relaxed to the new physics well 
before reaching the knee in the HR and log Toff— log g dia- 
grams. We check this by first doing the changes at different 
times in the post-VLTP evolution. We find that the tracks 
do not depend on the exact moment the changes are done, 
thus suggesting that the structure has already relaxed to 
the new situation. Additionaly, we estimate the thermal re- 
laxation time of the envelope as r ~ Sm*^^ CyTdm j — 
where Me.b. stands for the mass coordinate at the bottom 
of the envelope. We concentrate on the 0.53 and 0.584 Mq 
remnants of MA06. For these sequences we find that r is 
about one order of magnitude lower than the time it takes 
the remnants to evolve from log Tcff^ 4.6 to the knee in the 
HR diagram. In fact, r is about 2500 and 1600 yr for 0.53 
and 0.584 M© remnants respectively, as compared with the 
~ 23000 yr and ^ 12000 yr it takes the remnant to evolve 
to the knee. This guarantees that the envelope is thermally 
relaxed at that point. 



3.1. Equation of state and chemical composition of the C-0 
core 



To analyze the importance of the C-0 core composition and 
EoS we considered the 0.53 Mq post-VLTP sequence from 
MA06. With regard to the EoS we compare tracks result- 
ing from the use of the standard EoS of LPCODE — see 
Althaus et al. (2005) for references — with those coming 
from an updated version of Magni & Mazzitelli (1979) EoS. 
The latter is a more detailed equation of state which takes 
into account non-ideal corrections such as the pressure ef- 
fects on ionization and includes Coulomb interactions also 
in the non-degenerate regime — our standard EoS only in- 
cludes Coulomb corrections in the degenerate regime. To 
analyze the role of the core composition we reset the abun- 
dances below ln(l — m(r)/M*) — —2.4 — thus not altering 
the composition at the He-burning shell — to two extreme 
values: 92% of C and 92% of O by mass. This is not consis- 
tent with previous evolution but allows us to estimate its 
importance for post-VLTP tracks. 

Both changes in the EoS and the composition of the 
CO core do not yield significant changes in the post-VLTP 
tracks. In fact in none of these cases do we find the shift in 
log Teff to exceed 0.01 dex, being generally much smaller. 
Consequently, neither the C-0 core composition nor the 
EoS assumed in the computation of post-VLTP sequences 
play a role in the derivation of spectroscopical masses, and 
we can discard these two factors as possible reasons for a 
shift in post-VLTP tracks. 



3.2. Opacities 

Because the outer structure of PG1159 stars is completely 
ruled by radiative transport of energy, changes in the opac- 
ities could yield differences in the tracks. This may be par- 
ticularly interesting as the sequences of MA06 have been 
calculated for radiative opacities with a solar scaled metal- 
licity and PG1159 stars are known to present surface abun- 
dances rich in s-process elements and iron deficient (Miksa 
et al. 2002, WH06). In this regard, by how much the trans- 
formation of iron into heavier elements may alter the opac- 
ities in PG1159 is not known — for example, in a different 
context, Jeffery & Saio (2006) find differences in the pulsa- 
tion properties of subdwarf stars depending on whether it 
is iron or nickel that it is enhanced. Also the exact value of 
the original metallicity of the progenitor star of PG1159 is 
not known. We analyze the effect of changing both radia- 
tive and conductive opacities with very different results in 
each case. Full calculations of the VLTP and post-VLTP 
by means of consistent opacities are out of the scope of this 
work, however we can try to get an idea of how much the 
opacities affect the post-AGB tracks by artificially chang- 
ing the opacities in the post-VLTP evolution by arbitrary 
factors or by adopting different opacity tables. 

As a result of these experiments we find that for conduc- 
tive opacities even a change of 3 orders of magnitude do not 
produce significant changes in the post-VLTP tracks. Quite 
on the contrary, the tracks are more sensitive to radiative 
opacity changes. In fact we find that — for both the 0.53 
and the 0.584 Mq sequences — increasing the opacities by 
a factor 1.5 produces a reduction of ^ 0.04 dex in Tes and 
~ 0.2 dex in logL/L©. Similarly a reduction in the opacity 
by a factor 0.5 leads to increases of ~ 0.075 dex and ^ 0.3 
dex in temperature and luminosity, respectively. These are 
important changes and would clearly affect spectroscopic 
mass determination^. The shift in the location of post- 
VLTP tracks due to changes in the opacities is displayed in 
Table 2, where we show the change in log Tcff for different 
values of g and for two different remnant masses (0.53 and 
0.607 Mq). Also the induced shift in the mass derived from 
comparison with the g T^s values of MA06 tracks is shown. 
Two things deserves comments. The effect of different ra- 
diative opacities is much larger for higher remnant masses 
and at larger luminosities (i.e. lower gravities). Indeed, note 
that for the 0.53 Mq remnant an increase in the opacity of 
50% would not produce a shift of more than 0.01 Mq in 
spectroscopical mass determinations, and for the 0.607 Mq 
remnant the increase in the spectroscopical mass becomes 
very important, reaching up to 0.07 Mq at high luminosi- 
ties. Note also that the shift in log Toff is almost the same 
for the same change in k regardless of the mass. 

Due to the importance of this effect we consider inter- 
esting to analyze if the effect is due to the value of the 
opacity at some speciffic region of the star — e.g. the He- 
burning shell. We proceed then to make localised changes 
in the opacity and found, against expectations, that it is 
not the value of the opacity (per unit mass k) at some 
particular region that is relevant but the total opacity of 
the envelope (/g„veiopc i^dm). By looking at the models, we 

^ It is worth noting that we do not expect important changes 
in asteroseismological inferences due to changes in opacities. 
This is so because asteroseismological determinations are usu- 
ally based on adiabatic period studies, which are barely affected 
by changes in the opacities. 
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Table 2. Shifts in effective temperature ((51og Tcff ) induced by changes in n for different values of g. The value between 
brackets is the predicted induced shift in spectroscopical masses (in Mq). 



find that altering the radiative opacity produces almost no 
change in the structure of the envelope. Then, as dT / dm is 
not altered by changes in k, varying k leads to an opposite 
and proportional change in the luminosity Z(m) of the star 
via the relation 



dT_ 

dm 



kI 

,4y3 



(1) 



A clue of why only / reacts to a change in k can be obtained 
from the following simple analytical argument. If we assume 
that the envelopes of these objects are nearly homological 
to each other, then we have that under homology changes 
(with X = 5r/r) the change in the pressure and density of 
a shell is (see Kippenhahn & Weigert 1990 for a deduction) 



— = —4a;, — = —3a;. 
P p 



(2) 



Then if we assume an ideal gas equation of state for the 
envelope — which is quite correct — we have the additional 
relation 



5T_ 



P 



6P_ 



(3) 



Using the equation of the temperature profile (Eq[T|) by 
imposing an arbitrary change in the opacity 5k,/ n and using 



^. dT\ _ dT ^5T 
^dm J dm T 



we get 



5k 51 ^5r ^ 5T 
K I r T 



and by using Eq. 2 and 3 we finally find 

51 5k 
I k ' 



(4) 



(5) 



(6) 



which is quite similar to what it is observed in the numerical 
models — during the horizontal part of the tracks. Let us 
note that this change in the luminosity does at first order 
balance (in Eq.[T]) the change in opacity, leaving only second 
order effects on the factor kI to be balanced by the other 
factors in Eq. [1] 
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(7) 



Also, note that due to the high powers of r and T in Eq. 
[U small changes in these quantities should be enough to 
balance the remaining second order effects. In fact, when 



looking at the numerical models all r{m), P{m),T{m) re- 
main almost unchanged by the change in the opacity, being 
l(m) the only structure variable that undergoes an impor- 
tant variation. The change in l(m) seems to be associated 
with a change in the energy liberated by the helium burn- 
ing shell — change which can be attained with almost no 
change in T{m) due to the extremely high sensitivity of 
triple alpha reaction rates with temperature — and by a 
change in the dS/dt term of the energy equation — the 
lower the opacity, the faster the evolution and contraction 
of the envelope. Sumarizing we can say that altering the 
radiative opacity of the envelope leads to a similar change 
in the l{m) profile of the star which balances (at first or- 
der) the effects of the opacity change in Eq. 1. It seems that 
as consequence of this balance only minor changes appear 
in the other structure variables which remain almost un- 
changed — this probably reflects the fact that the run of 
these variables in the envelope of the star is forced by the 
radius and mass of the He- free core where most of the gravi- 
tational field is generated. Although this does not intend to 
be a complete explanation, something which is impossible 
nowadays due to the lack of an accepted explanation for the 
behaviour of structures with burning shells q we think that 
it sheds some light on what is happening on these models. 
Finally let us mention that as r(m) is not changed, then the 
radius of the star is almost the same, independently of the 
value of the opacity. Then, due to Steffan-Boltzmann's law 
we have that the change in the opacity produces a variation 
in the effective temperature of 



5Tcs 



cff 



15_L^ 
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4 



(8) 



In fact for the 0.584 Mq models with normal and enhanced 
(for a factor 1.5, ^ ~ 0.5) opacity we have that, in the 
knee, ^ ~ 0.34 and ^ - 0.09 and thus 4£i«- - • i 

' -ij Jeff -I off 

Note however that ^^ 



4 



0.34 is far from the value 0.5 in ^ 
(^ is close to 0.5 in the horizontal part of the track, i.e. 
log Toff< 4.8). 



3.3. Effect of different compositions on the opacities 

As models seem to be sensitive to the value of the radia- 
tive opacity in the envelope, we have analyzed how much 
opacities can change due to different adopted compositions. 

Firstly, we have assessed possible changes in the tracks 
due to changes in the total amount of metals in the models 



In fact the problem is in some way related with the long 
standing problem of why stars become red giants; see Faulkner 
(2005) for a recent discussion, and review, of this issue. 
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Fig. 5. Value of the opacity for different adopted composi- 
tions. 

— with the exception of C and O which are always kept con- 
sistent with the envelope abundance. We did this by using 
OPAL C- O- enhanced tables for Z = 0.01 and Z = 0.03 

— all of our previous sequences correspond to Z ^ 0.02. 
The change in metallicity was done well before evolution 
reached the PG1159 stage. We find that the resulting shifts 
in the tracks are barely noticeable. In fact, at the knee in 
the HR diagram the Z = 0.01 and Z = 0.03 tracks differ in 
effective temperature by only ~ 0.006 dex. 

Secondly, we have explored the use of opacities fully 
consistent with the abundances of the models. This is not 
a minor issue as Ne and N can be much larger than their 
solar scaled values — also Mg can reach values of 2%; see 
Werner et al. 2004. Specifically, we have used the tool at OP 
project website (Badnell et al. 2005) which allows to cal- 
culate opacities for arbitrary compositions. In this case we 
have not made any track calculation but instead we have 
just compared the opacities for a given model (i.e. for a 
given T and p profile). We compare first OP and OPAL 
opacities for the same imput composition. The result is 
shown in Fig. [5] Note that, for log T > 7, OP opacities 
tend to be about 5 to 7 % larger. This would probably in- 
troduce a small change of about 0.01 dex in the tracks. The 
inclusion of N and Ne — with abundances consistent with 
those displayed by the stellar models — markedly increases 
the opacity values below log T ~ 7.3, but almost no changes 
are present at higher temperatures where most of the mass 
of the envelope is stored, see FiglS] As PG1159 stars are 
supposed to be iron deficient due to s-process (WH06), we 
have analyzed the extreme case for which all the iron was 
changed into Ni. In this case, the opacity bump is located 
at larger T values, thus increasing opacity between log T= 
7.3 and 7.6. Although this change in the opacity is not 
enough to reconcile the discrepancy between spectroscopi- 
cal and asteroseimic masses of PG1159 stars, it is important 
to note that modifying the heavy metal distribution does 
introduce a change in the opacities at high temperatures. 
Indeed, as opacities increase with the atomic number of the 
elements — due to the increase in the possible atomic tran- 
sitions, Roger & Iglesias (1994) — it remains to be seen to 
what a degree an important increase in the content of very 
heavy metals due to s-process (both in the AGB and at 
the VLTP) can increase the opacity at the bottom of the 
envelope. Note that because of the higher ionization po- 



tentials, those elements are expected to affect opacities at 
much higher temperatures than do Fe or Ni. 

4. Other issues 

4.1. External consistency 

We have compared our tracks with the few models of mod- 
ern H-deficient post-AGB sequences available in the liter- 
ature, particularly the 0.604 Mq model of Herwig (2005) 
and the 0.632 Mq model of Lawlor & Mac Donald (2006). 




Fig. 6. H-deficient tracks of Herwig (2005) and Lawlor & 
Mac Donald (2006) compared with our post-VLTP tracks 
(with masses 0.87, 0.741, 0.664, 0.609, 0.585, 0.565, 0.542, 
0.530, and 0.515 Mq). Also a non-late helium flash (but 
H-deficient) track from O'Brien (2000) is shown for com- 
parison. 

Note from Fig. [6] that our models show a good agree- 
ment with the tracks of both authors. The agreement is 
remarkable despite the very different input physics and evo- 
lutionary history of progenitor stars considered by those au- 
thors. Indeed, Lawlor & Mac Donald (2006) models do not 
include any kind of overshooting prescription and Herwig's 
model is the result of an initially 2 Mq star model and thus 
with a distinct previous evolution than our 0.609 Mq se- 
quence which comes from a 3 Mq model. In addition, the 
EoS are different in all the cases. This supports the findings 
in Sects 2 and 3. For a quantitative inference, we estimate 
masses for those tracks by comparing their relative loca- 
tion with MA06 tracks. We derive masses of about 0.611 
Mq and 0.623 Mq for Herwig and Lawlor & Mac Donald 
sequences respectively — note that the resulting Herwig's 
track becomes bluer than ours, leading to slightly lower 
spectroscopical masses. In both cases the induced shift in 
spectroscopical masses would be lower than 0.01 Mq, thus 
reinforcing the robustness of the MA06 post-AGB tracks. 

4.2. Comparison with H-burning tracl<s 

As mentioned early in this work H-burners old tracks are 
cooler than MA06 tracks, thus leading to a much better 
agreement with asteroseismology (WH06). As we discussed 
in Sect 2, the difference in the location between old and new 
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tracks cannot be tracked back to a distinct previous evo- 
lution, with the exception of the VLTP event. It is worth 
noting that old models (Blocker 1995a, b) are based on the 
Cox & Stewart (1970) opacities, in contrast to new models 
that use OPAL and molecular opacities. In this connec- 
tion we feel important to recall that already Dreizler & 
Heber (1998) noted a shift of 0.03 M© between old Wood 
& Faulkner (1986) and O'Brien & Kawaler (as they appear 
in Dreizler & Heber 1998, here shown in Fig. 6) helium- 
burning tracks, where the latter make use of modern OPAL 
opacities. Appart from possible changes that could arise 
from the different opacities used in the calculations, the dif- 
ference in the tracks can be expected from the very fact that 
the Blocker's tracks are H-burners while VLTP are helium 
burners. In fact, we note from our sequences that there is 
a noticeable difference in H-burning-post-AGB and post- 
VLTP tracks for low remnant masses (^ 0.53 Mq). This 
can be seen in Fig [3 where we show our post-AGB H- 
burning tracks (with H-rich surface compositions) of 0.517 
and 0.53 Mq compared with the post-VLTP tracks of sim- 
ilar mass. Note that post-VLTP tracks are certainly bluer 
than their H-rich counterparts and that they are more com- 
pact in the high gravity region of the diagram (log g > 6.5), 
as found by MA06. Also note from the same figure that our 
0.517 Mq track is very similar to Blocker's 0.524 Mq track, 
strongly suggesting that the difference with Blocker (1995b) 
tracks at (very) low masses is because MA06 tracks are 
post-VLTP tracks. Using H-rich tracks to determine spec- 
troscopical masses for PG1159 stars will certainly influence 
the result, in particular the stellar masses for high-gravity 
PG1159 stars would be largely overestimated. 
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Fig. 7. Comparison between H-rich, H-burning tracks and 
H-deficicnt, Hc-burning tracks of low mass. It is clear from 
the figure that H-burners have lower Tcffand are much sim- 
ilar to Blocker's tracks. Using H-rich tracks to estimate 
PG1159 masses can lead to important overestimations in 
the low mass range. 



4.3. Comparison of post-LTP and post-VLTP tracks 

As mentioned early in this work, ^"^N-deficient PG1159 ob- 
jects are probably the descendents of LTP events. Thus, a 
priori one should be careful about using post-VLTP tracks 
for all PG1159 stars. In this context we now turn to an- 



alyze the question if there are systematic differences be- 
tween post-VLTP and post-LTP tracks. From figure 1 of 
Herwig 2001 it seems that there is no differece once the 
star enters the PG1159 stage. However we will now ana- 
lyze a wider range of masses. In Fig. 8 PG1159-tracks com- 
ing from VLTP and LTP are compared for similar rem- 
nant masses. In the upper panel LTP tracks with different 
H-abundances are compared with VLTP tracks of similar 
mass. The ~ 0.515 Mq tracks correspond to the sequence 
analysed in Althaus et al. (2007). In these sequences two dif- 
ferent post-LTP evolutions have been considered. The first 
in which the final surface H-abundance is norma|f| and a 
second in which due to mass loss episodes the whole H-rich 
envelope was eroded, exposing the He-rich intershell. Due to 
the absence of the H-burning shell in the second case it lies 
very close to the post VLTP tracks. The second experiment 
is also shown in the upper panel of Fig. 8 corresponds to 
an LTP sequence (0.543 Mq) in which the total H-content 
of the star was artificially diluted to different depths, thus 
leading to different final surface H-abundances. As can be 
seen once the star reaches the PG1159 stage, the lower the 
surface H-abundance the closer the LTP-tracks gets to the 
VLTP track of similar mass. Finally, in the lower panel of 
Fig. 8 post-LTP tracks of H abundances close to the usual 
detection limit are compared with VLTP tracks of similar 
mass. From that plot is clear that for surface gravities above 
log 5 = 6, where almost all PG1159 stars lie, VLTP tracks 
and LTP tracks with low H-abundances are similar. Then 
using post-VLTP tracks for spectroscopic mass determina- 
tions of LTP objects with no detectable H does not seem 
to introduce any systematic effect on the mass determina- 
tion. On the contrary using post-VLTP tracks for hybrid 
PG1159 stars may produce an important underestimation 
of the mass. 

5. Discussion and final remarks 

In the present work we have analyzed how uncertainties 
in the modeling of H-deficient post-VLTP remnants could 
afi^ect spectroscopic mass determinations of PG1159 type 
stars. In Sect. 2, inspired by a suggestion in WH06 we have 
analyzed the importance of previous evolution. As the cal- 
culation of each full sequence is extremely time consum- 
ing (both computational and human, as at some stages the 
models need hand interaction to converge them) we had 
to restrict ourselves to a limited region of the parameter 
space. Even in this case some conclusions can be drawn. 
Third dredge up alone does not seem to change the theoret- 
ical log TcS-logg-M locus and consequently its uncertainties 
can not affect spectroscopic mass determinations. On the 
other hand differences in mass loss rates (i.e. TP-AGB life- 
times) alter the location of the tracks, but only slightly. In 
fact our simulations show that even a reduction by a factor 
of 3 of TP-AGB lifetimes would not increase spectroscopic 
mass determinations by more than ^ 0.01 5Mq. Going even 
further we have shown that even extreme mass losses that 
produce low mass remnants {'^ 0.515 Mq) from very dif- 
ferent progenitors than those in MA06, does not introduce 
important shifts in spectroscopic mass determinations, be- 
ing only ~ 0.01 Mq. In an even more extreme limiting case 
we have computed a sequence, in which no overshooting was 

^ In this sequence, due to the very low remnant mass, the low 
intensity of the He-flash does not lead to any 3DUP. 
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Fig. 8. Comparison between post-LTP and post- VLTP 
tracks. Upper panel: Comparison between VLTP tracks 
and post-LTP tracks of similar mass but different H 
abundances. Lower panel: Comparison between post- VLTP 
tracks with post-LTP tracks that display surface hydrogen 
abundances, close to the detection limit. 



considered in the whole evolution, and found a very sim- 
ilar post- VLTP track than MA06. All these experiments 
suggest, contrary to the argument in WH06, that previous 
evolution only plays a secondary (and not very important) 
role in determining the theoretical log Tcn-log.g-Af locus. 
In light of these results we have discussed why the argu- 
ment presented in WH06 does not apply to post- VLTP se- 
quences. We showed that it is the HeFC (if any), and not 
the HFC, mass and radius that is important for post- VLTP 
tracks. In particular we find that the HcFC converges faster 
than the HFC in the mass-radius diagram. However, shell 
homology relations (as those used to derive the luminosity- 
mass-radius relation) should not be taken too seriously in 
these models, as they neglect the importance of the enve- 
lope and only relate the luminosity of the burning shells 
to the properties of matter in the burning shells and to the 
values of mass and radius of the core, since we find (in Sect. 
3.2) an important dependence of the shell luminosity with 
the whole opacity of the envelope. 

We have roughly addressed the robustness of the tracks 
regarding EoS, C-0 core composition, conductive and 
radiative opacities. We find that only radiative opacity 
may affect the location of the tracks to some an extent. 
Specifycally we find that the luminosity of the post- VLTP 



sequences in the horizontal part of the HR diagram is very 
sensitive to the envelope opacity. In fact the luminosity of 
the He-burning shell turns to be sensitive to the total opac- 
ity of the envelope. Wc also present some analytical argu- 
ments to explain the shift induced by changes in radiative 
opacitites. In this connection we explore how important the 
envelope composition can be for the opacity of the envelope. 
We find that changes in light metals (Ne and Mg) can make 
important changes in the opacities but only at low temper- 
atures (T < 5 X lO^K) where no much mass of the envelope 
is stored. Although this may be important for pulsational 
studies of PG1159 stars, it will certainly not change the 
Toff of the sequences. By contrast changing Fe into Ni in 
the opacity calculations we find a more slightly important 
change. This particularly leaves open the question of how 
much opacities at the bottom of the envelope can change if 
important amounts of Fe are transformed into very heavy 
metals by s-processcs. Wc can conclude that, unless there 
arc important changes in the abundances of very heavy el- 
ements due to s-process, an increase in the opacity at high 
T is not expected to change more than 10%. 

All these arguments show that MA06 tracks are robust 
enough as to be used for spectroscopical mass determina- 
tions of PG1159-type stars (specially at high gravities; log 
5^6). This robustness is reinforced by the good agree- 
ment (which corresponds to differences of <^ O.OIMq in 
spectroscopic mass determinations) between those tracks 
with the other modern post VLTP tracks available in the 
literature (Herwig 2005 and Lawlor & Mac Donald 2006). 
We have also addressed in Sect. 4.3 if any systematic in the 
mass determination may be due to the fact of some PG159 
stars being post-LTP objects with H-abundances below the 
detection limit. We find that the resulting tracks in the 
PG1159 region of the T^s-g (log g > 6) diagram are very 
similar to post- VLTP tracks when surface H abundance is 
below ^ 5% by mass fraction. Thus, wc conclude that the 
post- VLTP tracks of MA06 are solid enough for spectro- 
scopic mass determinations of post-LTP objects with H- 
abundances below the detection limit and, thus, it seems 
that no systematic should be present due to this effect. 
On the contrary, we find that using post- VLTP tracks for 
PG1159 stars with important H-abundances (the so called 
hybrid PG1159 stars) may lead to an important underesti- 
mation of the mass. Regarding the difference with Blocker's 
H-rich post-AGB tracks we can say that, for low mass rem- 
nants {<j O.53Af0), the differences in the tracks seem to be 
mainly due to the fact that those tracks are H-burners since 
our own H-burncr sequences are much colder than our post- 
VLTP ones. Other differences with older tracks may be due 
to the difference in the opacities adopted for the He, C - 
rich intershell (note that older works make use of old Cox 
& Stewart opacities). These seems to be supported by the 
good agreement between all H-deficient tracks that include 
modern OPAL opacities — Herwig 2005, Lawlor & Mac 
Donald 2006 and, more roughly (~ O.O2M0), even with 
the non- late helium flash 0.573 Mq sequence of O'Brien 
2000. 

From the present work we judge that the systematic 
discrepancy between asteroseismological and spectroscopi- 
cal mass determination methods should not be attributed 
to uncertainties in post-AGB tracks. Whether the discrep- 
ancy comes from errors in asteroseismological or spectro- 
scopical determinations is not known, however some points 
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are worth emphasising. Although asteroseismology is usu- 
ally accepted as a more accurate method (very low error 
bars are usually given), its robustness is not so clear. In 
fact recent works (Corsico & Althaus 2006, Corsico et al. 
2007a and Corsico et al. 2007b) show the results of aster- 
oseismology to be method dependent. In this context it is 
worth emphasising that the asteroseismic mass of PG 1 ISO- 
OSS is reduced to OMMq —only ~ 0.02 Mq higher than 
its spectroscopical mass — when detailed evolutionary mod- 
els and averanged period spacing (instead of the usually 
adopted asymptotic period spacing) are used in the analy- 
sis, see Corsico et al. (2006). Interestingly enough, during 
the referee stage of this article a new study of PG 0122-1-200 
(Corsico et al. 2007b) which is based on our evolutionary 
models and a detailed period by period fitting procedure, 
reduces the mass discrepancy (with MA06 value) in this 
star to less than a 4%. This clearly shows the existence of 
serious systematics in standard (i.e. based on asymptototic 
period spacing) asteroseismological determinations. In this 
context is worth noting that a mean PG1159 mass of 0.S73 
Mq like the one deduced from MA06 tracks, even if sensi- 
tively lower (0.044 Mq) than previously thought, is in good 
agreement with that of their probable descendants ^the DB 
white dwarfs (0.585 Mq, Beauchamp et al. 1996]3. Then 
our results not only call for a revision of PG11S9 model 
atmospheres but, specially, for a revision of systematics 
in usually adopted asteroseismological mass determination 
methods. 

Our full set of evolutionary tracks for post- 
VLTP objects is available at our web site at 
: / / www . f caglp . unlp . edu . ax/ evolgroup/. 
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